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ntEE VIBRATIONS OF THE ERDA-NASA 100 KW WIND TURBINE 


by C. C. Chamls* and T. L. Sullivan** 

Lewis Research Center 

ABSTRACT 

The ERDA-NASA wind tuitino (windmill), which consists of a 93-foot truss 
tower, a bed plate that supports mechanical and electrical equipment, and two 
02 5- foot long blades, was analyzed to determine its free vibrations using 
NASTRAN. The finite element representation of the system consisted of beam 
and plate elements . 'I'he free vibrations of the tower alone, the blades alone, 
and the complete system were determined experimentally in the field. These 
results were obtained by instrumenting the tower or blades with an accelerometer 
and impacting the components with an instrumented mass. 'Fhe predicted results 
for natural frequencies and mode shapes were in excellent agreement with meas- 
ured data. 
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FREE VIBItATIONS OF THE ERDA-NASA 100 KW WIND TURBINE 


by C. C. Chamis* and T. L. Sullivan** * 

Lewis Research Center 

INTRODUCTION 

In oixler to study some of the technical and economic problems involved in i 

generating electricity from wind energy, a 100 kW wind turbine generator (wind- 
mill) has been constructed at NASA-Lewis Research Center's Plum Brook 
facility near Sandusky, Ohio. The project is a cooperative effort between NASA 
and the Energy Research and Development Agency (ERDA). The wind turbine 
consists of a 93-foot tall square open truss steel tower; a bed plate that sits on 
top of the tower; generator, mechanical equipment and drive shafts which rest 
on the hpd plate; and twe 62. 5-foot long aluminum blades. 'Hie blades are 
downwind of the tower and inclined at an angle cf 7^ from the vertical plane. 

The wind tuiirine is designed to produce 100 kW of electrical power in an 
18-mph wind. 

Spera (ref. 3) performc“d a preliminary analysis comparing dynamic loads 
and stresses in two types of rotors for this wind turbine, under static, rated, 
and overload conditions. Blade vibrations were limited to the first beamwise 
mode. Closed-form equations were derived for the natural frequency in this 
mode. Donham, Schmidt, and Linscott present a detailed analysis and meas- 
ured data for the wind tuibine blade in reference 1. In the presen; investiga- 
tion the complete wind tuibine as well as individual components were analyzeii 
using NASTRAN to detennine .heir free vibrations. The free vibrations of the 
tower alone, one blade alone, the blades mounted on the bed plate, and the com- 
plete system were determined experimentally in the field by the Mobile Vibra- 
tion Lalxratoi*>’ of the tX'partment of Mechanical F'nginecring, Universit>' of 
Cincinnati (ref. 2), Briefly, these results were obtained by instrumenting the 
tower or blades with an accelerometer and impacting the components with an 
instrumented mass. Fast Fourier transform techniques were used to reduce 
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the measured data to free vibration results. 'ITie objective of this investigation 
v'as to verify the NASTIIAN finite element representation with measured data. 
Then the model can be used with confidence to fui*ther evaluate the dynamic 
response of the wind turbine. 

The finite element representation, the results obtained and comparison 
of these results with measured data are described and discussed in detail in 
this paper. 


DKSCRIPTION OF THE WIND TURBINE 

For puiposes of detcribing the wind turbine (fig. 1) it will be divided into 
three major sub-asserabl’es. 'Ihese are (1) tower, (2) bed plate and machinery, 
and (.1) blades. This section describes in detail these sub-assemblies. 


Tower 

The tower (fig. 2(a)) is a 93-foot tall welded open truss constructed from 
steel pipe (the four vertical columns), angle, channel, and wide flange beam 
sections. It is square and ranges from 30 feet wide at the base to 0.7 feet wide 
at the top. ITie tower weight is approximately 44,000 pounds. It is Ixilted to a 
steel -re info reed concrete foundation. Inside the tower a staii’way runs from 
ground level to the top ofthe tower. The stairway weight is approximately 12,000 
pounds. On the outside of one side cf the tower two beams run from ground level 
to the top of the tower, 'fhese serve as tracks for a self-propelled elevator 
(spider) used to transport tools and equipment from the base to the top oi the 
tower. 

The transition between the tower and the bed plate is made with a 4-foot 
rib-stiffened ronical transition section, the lower part i f which is bolted to the 
top of the tower. A motor driven gear drive within this cone orients the bed 
plate and blades with respect te the wind. The yaw plane is at the 97-foot level. 


Bed Plate and Machinery 

schematic of the bed plate and machinciy is shown in figure 2(b). I'he 
bed plate is a welded box beam and is attached through a large diameter bearing 


to the top of the conical transition secti'n. The bed plate supports the elements 
necessarj’ for transforming rotai 7 encrg>’ into electrical energy, 'lliese eU*~ 
ments are a hub, drive shafts, bearing supports, gearbox, generator, and 
hydraulic power supply 'ITie bed plate and machineiy weigh approximately 
.TO, out) pounds. 

The hub contains the hydraulically actuilcd gear mechanism that controls 
the blade pitch angle. It is connected to a drive shaft which is supported by 
two bearings. The center line of the rotor shaft is at the 100-foot level. The 
shaft is connected to a 15.1 gearbox to increase the rotational speed from 
40 to 1800 ipni. The high-speed ouiput shaM from the gear box is connected 
to the generator with V -belts. 


blades 

'Ihc blades (fig 2(c)) are fabricated from aluminum except for the root 
fitting whicii is steel. Ilie blades are bolted to the hub. Blade fabrication 
was based on aircralt structural tci4inolog\- Briefly, each lilade consists of 
thick aluminum leading edge skins and a large number of formed ribs that 
support the thiniK'r trailing edge skins. Kach blade conforms to a prescribed 
aerodynamic shape which has a large amount of twist ami taper. Total twist 
is ;i4®. Blade chord length varies 4 5 feet at the root t«* 2 feet at the tip. 
Blade thickness varies from 1 . •> feet at the root to 2 inches at the lip Kach 
blade weighs approximately 2000 pounds A more detailed description of the 
blades can be found in reference 1. 


DKSnUP riON OK niK NASTILW M(JI)K1, 

The windmill was modeled using NASTIb\N bi*am (('BOD and CBAB) and 
plate (C’(’rAD2) elements A schematic of the model is shown in figure .‘I. 

It consists of rtOT elements ami 170 notles with each no<le having 0 degrees of 
freedom. This section describes in detail how the system was moilcled. 


J 
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Tower 

The tower was modeled entirely with CROD and CBAR elements. The 
tower members were all common structural members and the section proper- 
ties used (moments of inertia about principal axes, torsional constant, and 
cross sectional area) were handbook values, 'llie tower model included the 
spider tracks that ran from the base to the top of the tower. The model in- 
cluded the weight but not the stiffness contribution of the stairw-ay. 


Top of Tower 

Plate elements (CQUAD2) in combination with bar elements (CBAR) were 
used to model the top of the tower. The thickness used for the plate elements 
was that of the actual structure. 


Conical Transition 

This section connects the top of the tower with the bed plate. It was 
modeled with plate elements (CGL’AD2). The element thickness was selected 
to provide bending stiffness equivalent to that cf the actual structure. 


Bcil Plate, Bearing Supports, and Drive Shaft 

The bed plate, bearing supports am. drive shaft were modeled with bar 
elements (CBAR). The properties for these elements were obtained by cal- 
culating the moments of inertia and cross sectional area for the various com- 
ponent sections. The weight of the machinery on the berl plate was included 
as nonstruciural mass in the appiopriate finite element. 


Blades 

The blades were modeled with bar elements (CBAR). Blade stiffness, 
twist and weight distribution are shown in figure 4. They were provided by 
the manufacturer of the blades 
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DISCUSSION OF RESULTS AND COMPARISON WITH MEASURED DATA 

NASTRAN vibration results for ore blade, blades mounted on bed plate, 
tower alone and the combined system are presented below and compared with 
measured data. O 


Blade 

The vibration free frequencies for one blade are summarized in Table I 
for the first five natural frtjquencies. Note that each frequency is identified 
as beamwise mode, chordwlse mode, or combined mode. The beamwise 
mode .shape is a predominantly o'lt-of-plane vibration motion with the plane 
taken parallel to the chord at the 3/4 blade span radius. "Chordwise" de- 
notes the direction parallel to this plane (li^. 4). The measured data shown 
in Table I were obtained by Lockheed-Califomia Co. and are reported in 
reference 2. Also included in Table I are analytical results from references 
1 and 3. As can be obsei*ved from the results in Table I, tlie NASTRAN 
predicted first natural frequency is within 2 percent of the measured value 
with "no” tip weight. 'Fhe corresponding higher frequencies (2 through 5) 
are within .5 percent. H is interesting to note that the predicted fourth and 
fifth natural frequencies are practically equal to the measured values with 
tip weight. This may be due to the fact that the blade has relatively large 
local variations along the span in both mass and stiffness that are not 
accounted for in the model. These large local variations tend to influence 
the higher frequencies especially when the wave lengths approach the length 
of the local variation. 

The predicted mode shapes for the first and second natural frequencies 
(first beamwise and first chordwise modes) are plotted in figures 5(a) and (b) 
w'here measured data (ref. 2) are also shewn. As can be obsened, the pre- 
dicted dynamic node shapes provide a good fit to the measured data. The 
predicted mode shape for the thiid natural frec(uency (second flap) is plotted 
in figure 5(c) where the ror responding measurtxi data are also plotted. For 
this case the predicted mode shape is in excellent agreement with the meas- 
ured data. 

The importarl conclusion to be drawn from the above discussion is that 
the finite clement ripresentation of the blade properly accounted for the 


6 


physical characteristics of the blade. This is a necessary requirement before 
the combined system can be modeled satisfactorily using finite elements. 


Blades Mounted on Bod Plate 

Before the blaoes and betl plate were lifted to the top of the tower, meas- 
ured data were obtained with the blades horizontal and in the powered posi- 
tion (ref. 2). These data are compared to NASTRAN result . for the same 
condition in Table II. F’or comparison purposes NASTIbXN results for the 
blades and hub only are includ»>d in this table. 

Characterization of the NASTRAN results as "second chcrdwise" and 
"third bcaniwise" is only approximate as these are again combined modes. 
Note that NASTRAN predicts two first beamwise (and subsequent) modes. 

In one, the blade tips more in opposite directions (cyclic) and in the other, 
the blade tips move in the same direction (ccllective) . There were insuffi- 
cient measured data to determine if a mode was cyclic or collective. 

The predicted and measured results are In g<y>d agreement. There were 
insufficient measured data to characterize the 12.0 Hz mode. The NASTRAN 
results indicated that this mode contained significant bed plate bending. 'The 
effect of the bed plate mass and stiffness en blade frequency can be seen by 
comparing the two sets of NASTitAN results, it can also be s<‘on that the 
bed plate effect te-nds t/i induce a slight reduction in the blade Irequencies. 


Towe r 

The vibration free Inquemcies cf the tower predicted by NASTRAN are 
summarized in Talile Hi where measured lield data (ref. 2) are also shown. 
'Fhe (N-S) and (K-VV) refer to tower directions as noted in figure 2(a). As can 
be observed fi*om the results in Table III. the predicted and meusui-ed data are 
in excellent agivement. 

The tower mode shaj-c cer responding to the first (K-\\’) natural frequency 
is shown in figure fi where measured daUi are also shown As can be observed 
from the plots in figure <>(b), the tower undergoes slight N-S motion as well. 
Note that the tower does n't t deflect latenilly (fig 6(a)) like a simple canti- 
lever. It deflects slightly in the reverse direction initially and then crosses 
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over at about elcva’i- n ?! feet. The jmporUint observation fron. the results 
In figure 6 is that the predicted mode shape is identical to that measured. 

The tower mode sh^pe for the second (E-W) natural frequency is plotted 
in figure 7. As can be obser/ed in figure 7(b) the predicted results show con- 
traction of the tower top in both directions. This contraction is the sum of 
both "rigid body" tiltmg of the tower top plane due t^ tower vertical motion 
and elastic deformation in the members connecting the tower legs it the top. 
The experimental data shows motion fnly along the E-VV direction. Th’® 
difference may be caused in part by (1) the stiffness of the stairway which 
is not accounted for in the NASTKAN mode] and which will tend tc support the 
stairs weight and (2) the lack of measure-d data lor tower vertical m Hion. The 
important point to be noted from figure 7 is that the predicted tower second 
(E-U^ mode shape is in reas« nably go< d agree ment with the field measured 
data. 

The important conclusion ir m the ab< ve discussion is that the NASTIIAN 
finite element represenUition of the tower pr , eriy accr>unts for the physical 
characteristics of the tower. Iliis again is a necessaiy requirement forol»- 
taining satisfactory results from the finite element representation of the 
combined system. 


Combined System 

'Phe NASTHAN vibration frequency results fur the c ombined system and 
those measured in the field (ref. 2) are summariz<Hl in Table IV for the blades 
in the horizontal and feathered position The orientati m ef the blades and 
shaft with respect to the towvr for thf results is shewm in figure M. .As can 
be ubsciwed fmm the results in this table, lor the m''st t.ar», the agreement 
between predicted and measured results is good. The lowest frequency pretlic- 
ted was 1. 59 liz for the blades h iizontal ami in thi p‘j\M>red position. 

There is some difference betsveen the pr*-dieted and measured mode shapes 
for the modes with frequencies bet" een 2 and :t llz The N.ASTIIAN mode 
shapes were primanly a e rnbinati -n f t aer fn-nding .uid blade chordwise 
motion. The first tw'o ncasured mode shaj es re first bending modes along 
the two diagi.nals. while the third mode was chaiMcterizid as the blade first 
chcrd.rise mode The diag*'nal motion is shown in figure 'i(O) f 'r the 2.2 llz 
mode It is compared to the 2 2S llz predict! d mode shape. It is possibly the 
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influence of the stairway stiffness that causes the difference between the pre- 
dicted and measured mode shapes although gootl agreement for the tower alone 
first bending mode shapes was obtainetl. Another possibility is an error in the 
mass distribution of the machinery on the bed plate Note in figure 8(a) the 
good agreeme. 1 t between predicted and measi:re<i mode shapes as observed 
from a south elevation. 

The important concii sion to be drawn from the discussion of the com- 
bined system results is that both structural analysis and dynamic testing 
methodologies have advanced to the point where complex structural systems 
can be analyzed and field tested readily, reliably, and economically. 


Sl’MMAHY AND CONCLI SIONS 

The predicted frequency and mode shape results obtained using NASTRAN 
for one blade alone, the blades niOuntcHl on the bed plate, the tower and the 
combined system wero in good agreement with the measured data. The finite 
element representation properly accounts for the combincKl .system's physical 
characteristics and can, therefore, be used with confidence to further evaluate 
its dynamic response. 'Phe results of this investigation demonstrate the versa- 
tility of NASTIIAN for deteiTnining the free vibrations of complex stroctures. 

In addition, dynamic testing methodologies for detennining free vibrations 
have advanced to the point where complex structural systesm can be field 
tested readily, reliably and economically. 
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TABLE I. - SINGLE BLADE FKEE FHECUENCIES, HZ 


NASTHAN 

Analysis (ref. .1) 

Analysis with no tip 
weight (ref. 1) 

Measured with tip 
weight (ref. 1) 

Measured with no 
tip weight (ref. 1) 


(a) Combined mode. 

(b) Beamwise mode. 

(c) Chordwise mode, 


1 

2 

1.7E(b) 

2.73(c) 

1.77(b) 


..65(b) 

2.3S(e) 

1.65(b) 

55(c) 

1.75(b) 

2. 66(c) 



TAB! 1 ‘. II. - FHEE FRECUENCIES OF BLADES MOUNTED ON BED 
PLATE, HZ (BI^DES HORIZONTAL IN POWERED POSITION) 


NASTRAn'^^ NASTRAN Measured 


Mode 

1st Beamwise | 

^ Cyclic 
[^Collective 

1st Chordwise | 

^ Cyclic 
[^Collective 

2nd Beamwise | 

^ Cyclic 
Collective 

2nd Chordwisej 

^Cyclic 

Collective 

3rd Beamwise 1 

^ Cyclic 
Collective 

Bed plate bonding 



(1 ) Blades anri hub onl\ 
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TABLK III. - TOWEH ALONE FHEE 


FREQUENCIES, HZ 



NASTRAN 

Measured^^^ 

1st Bending (N-S) 

■LTf) 

4.7 

1st Bending (E-W) 

5.19 

5. 1 

2nd Bending (E-V\') 

9.16 

9.4 

1st Torsion 

10.1 

10.5 


(1) Kef. 2. 


/ 


TABLE IV. - COMBINED SYSTEM FREE FRECU?:NCIES, HZ 
(B1J\DES HOIUZONTAL AND FEATHERED) 


NASTRAN 

Measured^^^ 

Mode 


Fre- 

Mode 

Fre- 



quency 


quency 

Blade 1st beamwlse (cyclic) 


1 . 69^^^ 






Blade 1st beamwise 

1.73 

Blade 1st beamwise (o'llective) 

1.7.3 




r 

2 15 

Tower diagoniil bending (SE-NW") 

2.1 

Combined tower 


2.28 

'Power diagonal bending (SE-NE) 

2 2 

HondiriK and blade 





Chordwlsc modes 


2 7.3 






Blade Is* chortlwise 

3.0 


s 

2.87 



Tower 1st torsion 


9 56 

'Power 1st torsion 

9.8 


(1) Ref. 2. 

(2) The lowest frequeney obtained for the combined system was 1. .59 Hz for the 

blades horizontal in p-'wered position 
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Figure 4. Wind turbine blade lotal stillness and weight distrtoulion. 
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